O 
O 

(N 

X> 

IX. 
(N 



> 

in 

(N 
O 
<N 
O 
VO 
O 

22 

Oh! 

Oh! 
CD 
43 



X 



BU-HEPP-06-04 
Jan., 2006 



New Applications of Resummation in Non-Abelian Gauge Theories: 
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We present the elements of three applications of resummation methods in non-Abelian gauge theories: (1), 
QED®QCD exponentiation and shower/ME matching for LHC physics; (2), IR improvement of DGLAP theory; 
(3), resummed quantum gravity and the final state of Hawking radiation. In all cases, the extension of the YFS 
approach, originally introduced for Abelian gauge theory, to non-Abelian gauge theories, QCD and quantum 
general relativity, leads to new results and solutions which we briefly summarize. 



1. Introduction 

Multiple gluon (n(G)) effects are already needed in 
many high energy collider physics scenarios, such as tt 
production at FNAL, polarized pp processes at RHIC, 
bb production at FNAL,- ■ ■ , and contribute a large 
part of the current uncertainty on nit ,~ 2-3 GeV [1]. 
For the LHC, and any TeV scale linear collider, the 
rather more demanding requirements make exponen- 
tiated soft n(G) results, exact through Q(pt%), realized 
by MC methods in the presence of parton showers 
without double counting on an event-by-event basis 
with exact phase space, an essential part of the nec- 
essary theory. 

For such precision QCD results, ~ 1%, may QED 
higher order corrections also be relevant? The stan- 
dard treatments of resummation in the two theories 
treat them separately, where we have considerable lit- 
erature on resummation in QCD [2-5] and on expo- 
nentiation in QED [6-9]. We use the words resum- 
mation and exponentiation interchangeably here be- 
cause in all cases a leading exponential factor is iden- 
tified in the resummations which we discuss in this 
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paper. Results from Refs. [10-14] show a few per 
mille level QED effect from structure function evo- 
lution when QED and QCD DGLAP [15] kernels are 
treated simultaneously. After reviewing the YFS the- 
ory and our extension of it to QCD in the next Sec- 
tion, in Section[2]we combine [16] the respective two 
exponential algebras at the level of the amplitudes 
to treat QCD and QED exponentiation simultane- 
ously, QED$§QCD exponentiation, to be able to iso- 
late, on an event-by-event basis, the possible interplay 
between the large effects from soft gluons and soft 
photons. We also show that the QCD exponentiation 
extends naturally to quantum gravity and allows us 
to develop UV finite resummed quantum gravity [17]. 
We find surprisingly that our extension of the ideas of 
Ref. [7] to Feynman's formulation [18,19] of Einstein's 
theory, the UV behavior of the respective theory be- 
comes convergent. 

The size of the QED(g)QCD threshold effects is il- 
lustrated in Sect. |1] wherein we present as well a new 
approach to shower/ME matching [20] for exact fixed- 
higher-order results. We also present IR-improved 
DGLAP theory [21] in this Section. In Sect. we 
discuss the final state of Hawking [22] radiation for 
an originally very massive black hole in resummed 
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quantum gravity. 

2. Review of YFS Theory and Its Exten- 
sion to QCD 

We consider first the QED case presented in 
Refs. [23] and realized by MC methods, for which, 
as we illustrate here for e + (pi)e~(gi) — > f(p 2 )f{q 2 ) + 
n(7)(fei, •, k n ), renormalization group improved YFS 
theory [24] gives 

n = ' J j = l i 
d 4 V iy(p 1 +q 1 -p 2 -q 2 -Y.ik i ) + D (1) 

(2^)4 

d 3 p 2 d 3 q 2 



Pn(kl 



Pa 9a 



where the YFS real and virtual infrared functions 
B, D, B are known [7,23]. Examples of the YFS 
hard photon residuals ft in Q, i = 0, 1,2, are given 
in the fourth paper in Refs. [23] for the MC BHLUMI 
4.04. 

In Ref. [4] we have extended this YFS theory to 
QCD: 



SUM IR (QCD) ' 



E/n^ 



n=0 j=l J 
= 1 y(Pl + P 2 -Qi-02-E fcj) + -D QCD 



d*y 
(2^) 4 



(2) 



■Pn(kl 



i kn) 



d 3 P 2 d 3 Q 2 

e 2 Q 2 



the general exponentiation of the eikonal cross sec- 
tions for non-Abelian gauge theory; his result is an 
approximate one in which everything that does not 
eikonalize and exponentiate is dropped whereas our 
result is exact. 

3. Extension to QED®QCD and Quantum 
Gravity 

Simultaneous exponentiation of QED and QCD 
higher order effects [16] gives 



B 



jjnls I r^nls 

E>OC.r> + DflBD = -t>i 



QCD ~ > D QCD T DQED 



QCED, 



^QCD 



which leads to 



^QCD T D QED = OQCED, 



(3) 



00 /* n ,3 , m ,3 , / 
o SUM IR (QCED) ^ Y[ II J ' 2 



/ 



ri,m — 31=1 ^ 1 



32 = 1 



d y iy(pi+qi-y 2 -<l2-Y. t 'i 2 ) + I) QCED 

(2*)< 

J ffa k -ti k' ) d3p2 d3<12 

p 2 g 2 



(4) 



where the new YFS residuals /3 n ,m(ki, . . . , fc n ; fc^, . . . , k' m ), 
with n hard gluons and m hard photons, represent 
the successive application of the YFS expansion first 
for QCD and subsequently for QED. 
The infrared functions are now 

SUMm(QCED) = 2a s XB$S ED + 2a a B n Q u CED 



-Dqced = Dqcd + D. 



(5) 



where now the hard gluon residuals j3 n (ki, . . . , kn) 
represented by 



@n (kii ■ ■ - ; k n ) — ^ ^ P n (ki 7 . . . , k n ) 



are free of all infrared divergences to 
all orders in ct s (Q). The functions 

SUM ra (QCD), Dqcd, f3 n {kx, . . . ,k n ) are the 
QCD analoga of the QED functions 2a Re B + 
2a B, D, p n (ki, . . . , kn) and are defined in Ref. [4,5]. 
For reference, the respective process exponentiated 
in J5J can be taken to be q + q — > q" + q'" + n(G) 
for quarks q, q" and anti-quarks q',q"'. 

Our exponential factor corresponds to the N — 1 
term in the exponent in GatheralPs formula [25] for 



We have a new IR Algebra(QCED) which is char- 
acterized by the following estimates of the respective 
average values of Bjorken's variable: x a v g (QED) = 
j(QED)/(l+i(QED)), x avg (QCD) Si 7 (QCZ>)/(1+ 
y(QCD)), with 7 (A) = 2£aCa (Ls — 1), A = 
QED, QCD, Ca = Q 2 f ,C F , respectively, for A = 
QED, QCD, where L s is the big log, so that QCD 
dominant corrections happen an order of magnitude 
earlier than those for QED; the leading /3Q°Q°'-level 
gives a good estimate of the size of the effects we 
study. 

We can apply @ to quantum general relativity: 
For the scalar 2-point function we get [17] 



rcsummcd 



S' s + it) 



(6) 
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for 



2 4 f d 4 e 



1 



16tt 4 P 
1 



A 2 + it 



(7) 



[l 2 + 2£k + A + ie) 2 

This is the basic result. 

Note the following: £!, starts in 0(k 2 ), so we 
may drop it in calculating one-loop effects; explicit 
evaluation gives, for the deep UV regime, B'J( k ) — 

K g |^ I In ^ m ^| fc ^| J , s ° that @ falls faster than any 

power of \k 2 \; and, if m vanishes, using the usual — /i 2 



In 



normalization point we get B g (k) = * 

so that (jSJ again vanishes faster than any power of 
|fc 2 |! These observations mean [17] that, when the re- 
spective analoga of © are used, one-loop corrections 
are finite. Indeed, all quantum gravity loops are UV 
finite [17]. We refer to this approach to quantum gen- 
eral relativity as resummed quantum gravity (RQG) . 

4. QEDigiQCD Threshold Corrections, 
Shower/ME Matching and IR-Improved 
DGLAP Theory at the LHC 

We shall apply the new simultaneous QED(g>QCD 
exponentiation calculus to the single Z production 
with leptonic decay at the LHC ( and at FNAL) 
to focus on the ISR alone, for definiteness. See also 
the work of Refs. [26-28] for exact 0(a) results and 
Refs. [29-31] for exact 0(a 2 B ) results. 

For the basic formula 

da exp (pp -> V + X -> W + X') = 
^ / dxidx j Fi(x i )Fj(xj)da exp (x l Xjs), (g) 



we use the result in @ here with semi-analytical 
methods and structure functions from Ref. [32]. A 
MC realization will appear elsewhere [33]. 
Here we make the following observations. 

• In © we do not attempt at this time to replace 
Herwig [34] and/or Pythia [35] - we intend to 
combine our exact YFS calculus, da exp (xiXjs), 
with Herwig and/or Pythia by using them/it 
"in lieu" of {Fi} as follows: A. Use a Her- 
wig/Pythia shower for pr < /ti, and YFS nG 
radiation for pr > 

B. Or, expand the Herwig/Pythia shower for- 
mula (SdtJexp and adjust /3 n ,m to exactness for 



the desired order with new /3^ m . 
In either A or B, we first use {Fi} to pick 
(#1,312); make an event with da exp ; then shower 
event using Herwig/Pythia via Les Houches 
recipe [36]. 

• This combination of theoretical constructs can 
be systematically improved with exact results 
order-by-order in a s ,a, with exact phase space. 

• Possible new parton showers such as one based 
on the recent alternative parton evolution algo- 
rithm in Refs.' [37] can also be used. 

• Due to its lack of color coherence [38] Isajet [39] 
is not considered here. 

With this said, we compute , with and without 
QED, the ratio r exp = cr exp /aBorn to get the results 
(We stress that we do not use the narrow resonance 
approximation here.) 



1.1901 , QCED = QCD+QED, LHC 
1.1872 ,QCD, LHC 

1.1911 , QCED = QCD+QED, Tevatron 
.1.1879 ,QCD, Tevatron. 

0) 



We note that QED is at .3% at both LHC and FNAL, 
this is stable under scale variations, we agree with 
the results in Refs. [26,27,29,30], and the QED effect 
similar in size to the structure function results [10- 
14]. DGLAP synthesization [40] has not compromised 
the normalization. 

With the precision tag needed for the LHC in mind, 
we apply QCD exponentiation theory to DGLAP ker- 
nels [21]: we get 



° qq {z) = CfFy 



1 — z 

-f q ( lq )S(l-z)] (10) 



where f q (-y q ) = 
0o ' ~~ 2 



_2 2_ 

7<i 7g + l 
I a s Cf / 7r 2 

"T 7T \ 3 



1 

|) and 



r(i+7,) 



7, = C F ^t = 

F Y Fs{lq) = 

Similar results hold for PGq, Pgg, PqG, so 



4 



B.F.L. Ward 



that we have: 



P Gq {z) = C F FYF S ('y q )e2 l 



(1 



(11) 



Pgg(z) = 2C G F Y Fs(~fG)ei SG {- 



-{1-zT 



1+7G 



(1-z) 



+ z(l-z 



, 1 +7G 1 



f G (lG)5(l-z)}, 



P qG (z) = F YFS (-i G )e^ G -{z\l - z) 

2 "IG 1 



+ (l-z)V G }, 



where 7g = C G —t 
and 



00 



■, <5g 



7G I a s^G ( TI 
2 "T" 7T V 3 



(12) 



(13) 
-I). 



/g(7c 



n/ 



+ 



C G (1 + 7g)(2 + 7g)(3 + 7g) 
2 



7g(1 + 7g)(2 + 7g) 
1 



(1 + 7g )(2 + 7g) 2(3 + 7g)(4 + 7g ) 

+ (2 + 7g )(3 + 7g)(4 4-7g)' (U) 

Applying this new kernel set to the evolution ol the 
parton distributions (Recall that moments of kernels 
logarithmic exponents for evolution.), we have for 
the NS anomalous dimension A„ the result 

^ s = C F F YFS (j q )e^ 

[B(n, lq )+B(n + 2,~f q )-f q (~, q )] (15) 

where B(x,y) is the beta function given by B(x, y) = 
T(x)T(y) /T(x + y). Compare the usual result 



C F [- 



l 2 ? I] 

n(n + l) f^j' 



(16) 



The IR-improved n-th moment goes for large n to 
a multiple of —fq, consistent with lim n ^oo z n ~ 1 = 
for < z < 1; the usual result diverges as 
—2C F In n. The two results differ for finite n as 
well: we get, for example, for ce s = .118, A^ S = 
Cf(-1.33),Cf(-0.966) for (Q| and (O, respec- 
tively. For the n-th moment of the NS parton dis- 
tribution itself, M" S , we have [21] 

M^ s (t) = M^ s (t )e^ lE ^ s ^ to) '>- E ^ s ^ t))] 

t,t large with t>>t \ Ots(t) 




fc + q 
(a) 




9 r^A. 

+ <7 
(c) 



Figure 1. The graviton((a),(b)) and its ghost ((c)) 
one-loop contributions to the graviton propaga- 
tor, q is the 4-momentum of the graviton. 

(17) 

where Ei(x) — dre r /r is the exponential inte- 
gral function, 

2Cf 2? 
a n = -^—F YF s(^ q )e 4 [-8(71,7,) 
Po 

+ B(n + 2, 7 ,)-/,(7,)] (18) 

- _ , Ji (a a (to) - a a (t)) 
a n -a n[ L+ 2 ln(as(to)/Qs(t)) ) 

with 5i — |V This is to be compared with 

the un-IR- improved result where last line in eg. 1171 
holds exactly with a'„ = 2A„ S /j3o- Comparison 
with the higher order exact results on the kernels in 
Refs. [41,42] is done in Ref. [21]. 

5. Final State of Hawking Radiation 

Consider the graviton propagator in the theory of 
gravity coupled to a massive scalar(Higgs) field, as 
studied by Feynman in Refs. [18, 19]. We have the 
graphs in Fig.^and [5] Using the resummed theory, 
we get that the Newton potential becomes [17] 



, , GnM . 
*jv(r) = (1 



for a = 0.210A/ P ,. Our results imply 



G(k) = G N /(1 + ^) 
or 



(19) 
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the lapse function is, from Ref. [44], 



k + q 
(a) 



+ 




(b) 



Figure 2. The scalar one- loop contribution to the 
graviton propagator, q is the 4-momentum of the 
graviton. 



so that we have fixed point behavior for k — » do, 
in agreement with the phenomenological asymptotic 
safety approach of Refs. [43-45]. We can also see [17] 
that our results imply that an elementary particle 
has no horizon which also agrees with the result of 
Ref. [44] that a black hole with a mass less than 
M cr ~ Mpi has no horizon. The basic physics is 
the following: G(k) vanishes as 1/fc 2 for k 2 —* oo. 

There is a further agreement: final state of Hawk- 
ing radiation of an originally very massive black hole 
because our value of the coefficient, \, of k 2 in the 
denominator of G(k) agrees with that found by Bon- 
nanno and Reuter(B-R) [44], if we use their prescrip- 
tion for the relationship between k and r in the regime 
where the lapse function vanishes, we get the same 
Hawking radiation phenomenology as they do: the 
black hole evaporates in the B-R analysis until it 
reaches a mass M cr ~ Mpi at which the Bekenstein- 
Hawking temperature vanishes, leaving a Planck scale 
remnant. 

What is the remnant's fate? In Ref. [46] we show 
that our quantum loop effects combined with the G(r) 
of Bonnanno and Reuter imply that the horizon of the 
Planck scale remnant is obviated - consistent with 
recent results of Hawking [47]. To wit, in the metric 
class 



_ 2G(r)M 
r 

B{x) 
B(x) + 2x 2lx 



(21) 



: G N M ' 



where B(x) 



pi 



x 3 - 2x 2 + Qx + -yQ, for Q = 

G -jjj = M a - . After Hawking- radiating to regime 
near M cr ~ Mpi, quantum loops allow us to replace 
G(r) with Gat(1 — e~ ar ) in the lapse function for 
r < r>, the outermost solution of 



G(r) = G N {l~e- ar ). 



(22) 



ds 



f(r)dt 2 - /(r)" 1 ^ 2 



r 2 dn 2 



(20) 



In this way, we see that the inner horizon moves to 
negative r and the outer horizon moves to r = at the 
new critical mass ~ 2.38Mp;. We note that the re- 
sults in Ref. [48] show that loop quantum gravity [49] 
concurs with this general conclusion. 

We have arrived at the following prediction: there 
should energetic cosmic rays at E ~ Mpi due the 
decay of such a remnant. 
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